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ABSTRACT 

We present aperture synthesis imaging of dense molecular gas in the Large Magellanic Cloud, taken with 
the prototype millimeter receivers of the Australia Telescope Compact Array (ATCA). Our observations of 
the N113 H II region reveal a condensation with a size of ~6" (1.5 pc) FWHM, detected strongly in the 
1-0 lines of HCO + , HCN and HNC, and weakly in C2H. Comparison of the ATCA observations with single- 
dish maps from the Mopra Telescope and sensitive spectra from the Swedish-ESO Submillimetre Telescope 
indicates that the condensation is a massive clump of ~10 4 M Q within a larger ^10 5 M Q molecular cloud. 
The clump is centered adjacent to a compact, obscured H II region which is part of a linear structure of radio 
continuum sources extending across the molecular cloud. We suggest that the clump represents a possible site 
for triggered star formation. Examining the integrated line intensities as a function of interferometer baseline 
length, we find evidence for decreasing HCO + /HCN and HCN/HNC ratios on longer baselines. These trends 
are consistent with a significant component of the HCO + emission arising in an extended clump envelope and 
a lower HCN/HNC abundance ratio in dense cores. 
Subject headings: Magellanic Clouds — ISM: molecules — ISM: structure — galaxies: ISM 



1. introduction 

Cold molecular gas in the interstellar medium (ISM) is usu- 
ally studied via rotational transitions of trace molecules that 
emit at millimeter and sub-millimeter wavelengths. The most 
commonly used tracer is CO, which typically traces gas den- 
sities of ^10 2 -10 3 cm" 3 . Higher densities are best probed 
using molecules with larger dipole moments such as HCO + 
and HCN. Emission from these molecules comes from the 
densest condensations within molecular clouds, often referred 
to as "clumps" and "cores" (e.g., Blitz 1991). A molecu- 
lar core has typical radius R ~ 0.1 pc, mass M ~ 1-10 Af©, 
and density n > 10 4 cm" 3 , and is thought to eventually form 
a single star or multiple star system. Clumps are somewhat 
larger structures which may evolve into star clusters. The use 
of high-density tracers to study the initial conditions for star 
formation has long been common practice for Galactic ob- 
servers, and has recently gained greater prominence in stud- 
ies of nearby galaxies following the demonstration of a tight 
correlation between HCN and far-infrared (FIR) luminosities 
(Gao & Solomon 2004b), the latter tracing recent star forma- 
tion. 

Although they have similar dipole moments, HCO + and 
HCN may not trace dense gas equally well, as one must con- 
sider the chemical processes which lead to their formation and 
destruction. This becomes especially apparent when examin- 
ing molecular abundances in the Magellanic Clouds, where a 
lower dust-to-gas ratio and locally active star formation ac- 
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tivity combine to yield a stronger photoionizing flux than in 
the Galaxy. As a result, molecular clouds are immersed in 
extensive photon-dominated regions (PDRs) where carbon is 
mostly ionized, as evidenced by the much higher [C II]/CO in- 
tensity ratios in the LMC (Mochizuki et al. 1994; Israel etalJ 
1996). While the formation of HCO + in regions shielded from 
far-ultraviolet (FUV) radiation is dominated by reactions be- 
tween H3 and CO, it can also occ ur in PDRs via reacti ons be- 
tween C+ and OH, 2 , and H z O JGraedel et aljfl982l) . Thus, 
the spatial extent of the HCO + emission may be much larger 
than that of HCN — even comparable to that of CO. 

The HNC molecule, on the other hand, has been suggested 
to be a tracer of cold cores, as steady-state chemical mod- 
els predict that it will be much weaker than HCN in regions 
of high temperature and density (Schilke et al. 1992). Con- 
sistent with these predictions, the HCN/HNC abundance ra- 
tio is particularly high («80) in the hot core of Orion KL 
(Schilke et al. 1992), whereas i t is close to u nity in cool dark 
clouds (Churchwell et al. 1984J: lHirota et alll 998). However, 
observations of starburst galaxies, where temperatures are ex- 
pected to be high, have revealed HNC intensities th at are 
sometimes close to those of HCN ( Huettemeister et al. 1995; 
Aalto et al. 2002). While an extended component of cool, qui- 
escent gas could contribute to the HNC emission, Aal to et alJ 
(2002) speculate that PDR chemistry may play a role: re- 
actions involving HCNH + may become important, which re- 
combines to produce HCN and HNC with equal probability 
and largely independent of kinetic temperature. Thus, HNC 
is potentially another tracer of PDR conditions. 

As one of the initial targets for the prototype millimeter re- 



ceivers on the Australia Telescope Compact Array 2 (ATCA), 
we have observed the H II region Nl 13 in the Large Magel- 
lanic Cloud (LMC) in the 3mm transitions of HCO + , HCN, 
HNC and C 2 H. The LMC's proximity (d=52 kpc, Panagia 
1999) and low metallicity make it a very attractive target for 
studies of the molecular ISM, but its southern declination 
has heretofore made it inaccesible to millimeter interferom- 
etry. We use HCN as a fiducial dense gas tracer and examine 
the relative intensities of the other molecular lines (including 
C2H, produced in reactions involving C + ) for enhancements 
due to the presence of a PDR. In particular, we examine the 
line ratios as a function of spatial scale, which allows a com- 
parison of compact and extended emission regio ns. 

The molecular cloud associated with N113 (Henize 1956) 
is one of the richest in the LMC, with a peak CO(l-O) 
brightness temperature (measured with the 45" beam of 
the Sw edish-ESO Submillimetre Telescope, SEST) of —9 
K dChinetal.1119971) . Also associated with N113 are three 
young (<10 My r) stellar clusters NGC 1874, NGC 1876, 
and NGC 1877 JBicaet alj fl992). Current star formation 
activity appears concentrated in three compact radio contin- 
uum sources, superposed on an extended emission component 
and aligned in a northwest-southeast direction with separa- 
tions of —20" (Brooks & Whiteoa kTl997h . The faintest, east- 
ernmost source contains by far the brightest H2O m aser in 
the LMC (Whiteoak& Gardner 1986; Lazendic et al. 2002) 
as well as an OH maser ( Brooks & Whiteoak 1997). The 
bright est, westernmost so urce also shows H2O maser emis- 
sion {Lazendic et al. 2002J); it is towards this source that we 
have targeted our ATCA observations. All four of the 3mm 
transitions observed with ATCA have been previously de- 
tected in this region with the SEST. In this paper, we present 
the ATCA results and compare them with single-dish results 
from the Mopra and SEST telescopes. 

2. OBSERVATIONS 

2.1. ATCA 3-mm data 

Observations with the ATCA were made in 2003 July and 
August in the EW214 and EW367 configurations respectively. 
At the time of the observations, the ATCA had three anten- 
nas of 22 m diameter equipped with dual polarisation 3-mm 
receivers covering the bands 84.9-87.3 and 88.5-91.3 GHz. 
At these frequencies the primary beam has a FWHM of 36". 
Both configurations had all three antennas arranged in an east- 
west line. The EW214 configuration provided baselines of 
30m, 75m, and 105m, whereas the EW367 provided baselines 
of 45m, 90m, and 135m. The pointing center for the obser- 
vations was a 2 ooo=5 h 13 m 18 s , <5 20 oo=-69°22'20". The ATCA 
correlator can be configured to observe two frequencies si- 
multaneously, and we cycled between 3 pairs of frequencies 
during each observing session in order to obtain comparable 
coverage of the visibility plane. The first pair of frequencies 
consisted of the HCO + (7=1-0) and HNC (7=1-0) lines at 
89.1885 and 90.6635 GHz respectively, redshifted to the ap- 
propriate observatory frequency for vlsr=235 km s" 1 , the sec- 
ond pair consisted of the HCN (7=1-0) line at 88.6318 GHz 
and HNC, and the third pair consisted of the C 2 H (N=l-0, 
7=3/2-1/2, F=2-l) line at 87.3169 GHz and a 128 MHz wide 
continuum window at 86.243 GHz. Each line was observed 
with a channel spacing of 0.5 MHz across a 64 MHz band- 
width, except for C 2 H, which was observed with a channel 

2 The Australia Telescope is funded by the Commonwealth of Australia 
for operation as a National Facility managed by CSIRO. 
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TABLE 1 

Properties of the ATCA images 

Line Beam size Beam PA cr ch a <j cll T™,^ Vel. range Flux b 

(" X ") (°) (mJybnT 1 ) (K) (K) (kms' 1 ) (Jykms' 1 ) 

0.38 4.2 [229,241] 5.2-8.8 
0.39 2.8 [223,243] 4.1-6.9 
0.21 1.4 [229,241] 1.0-2.1 
0.46 1.1 [233,239] 0.5-1.5 
0.03 0.3 — 0.04 Jy 

a RMS noise for a 2 km s channel (line images) or continuum integrated 
across 84 MHz. 

b Flux measurements are made in boxes of 10" X 10" and 20" X 20" cen- 
tered on a2000=5:13:17.2, <52ooo= - 69:22:23. 



spacing of 0.0625 MHz across 16 MHz. Figure [Q a) shows 
the coverage of the visibility plane obtained for the HCO + 
line. 

The da ta were cal ibrated and imaged with the MIRIAD 
package (ISault et al.l 119951) . which at ATNF has been up- 
dated with new procedures to handle ATCA 3-mm data. Flux 
calibration was performed using Uranus, while the quasars 
B0537-441 and B 1253-055 (3c279) were used for gain and 
bandpass calibration respectively. However, for the last fre- 
quency pair (C 2 H and continuum), the SiO maser R Dor was 
used for gain calibration instead; in doing so the maser flux 
was assumed to be constant over the course of an —8 hour 
observation. An instrumental phase error causing phase dis- 
continuities when changing sources had been identified with 
ATCA 3-mm data, 3 and we compensated for this using an an- 
tenna position correction which was determined by observa- 
tions of multiple quasars within a few days of the observa- 
tions. We produced channel maps spaced in velocity by 2 
km s" 1 , which is close to the effective resolution provided by 
the correlator in the HCN, HNC, and HCO + lines. The phase 
center was shifted to a 20 oo=5 h 13 m 17?2, 5 200 o=-69 22'23" be- 
fore imaging, a shift of 5" (this shift has negligible impact 
on the primary beam response at the locations of detected 
emission). The noise level and synthesized beamwidth for 
each map are given in Table ^ After imaging the maps were 
CLEANed to a 2rr level. 

2.2. ATCA 1,3-cm data 

Additional observations were obtained with ATCA in 2005 
March using the recently installed 12-mm receivers. N113 
was observed in the H214 configuration (baseline lengths 
80m-250m) over two days with a total integration time of 
160 min. The correlator was configured for 512 channels 
across 8 MHz (velocity resolution 0.24km s" 1 ) at 23.673 GHz 
and 32 channels across 128 MHz at 23.9 GHz. The nar- 
rowband observations were designed to detect the NH3(1,1) 
transition but failed to do so at an RMS noise level of 20 
mJybm" 1 . The wideband observations, used to image the 
continuum, achieved a noise level of —0.4 mJy bm" 1 . Flux, 
bandpass, and phase calibration were performed using the 
quasars B 1934-638, B1921-293, and B0454-8 10 respectively, 
following standard reduction procedures in MIRIAD. After 
Fourier inversion the 1 .3 -cm continuum image was CLEANed 
and restored using a beam size of 10"1 x9."l. Note that be- 
cause of the limited coverage of the Fourier plane achieved 
with the H214 configuration, structure on scales much larger 

3 This is believed to be an instrumental error which corrupts the antenna 
baseline solution. 
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FIG. 1 . — Summary of ATCA observations towards Nl 13. (a) Coverage of the visibility plane for HCO + obtained from combining both runs. Other panels 
show integrated intensity maps for (b) HCO+; (c) C 2 H; (d) HNC; (e) HCN; (f) continuum. Offsets are relative to a 200 o=5 h 13 m 17?2, <5 2 ooo=-69°22'23". The 
synthesized beam FWHM is shown in the lower left corner and the large circle represents the FWHM field of view. The solid square and cross represent the 
positions of the A'-band source and the weaker H 2 maser respectively. The contour levels are spaced by 6 K km s _1 for HCO + and HCN, 3 K km s for HNC 
and C 2 H, and 0.06 K for the continuum. 



than 30" will be poorly imaged. Thus, the resulting image 
is dominated by compact ionized regions, and the extended 
continuum emission from the large-scale H II region was not 
detected. 

2.3. SEST data 

Although most of our comparison with SEST data m akes 
use of values previously presented by Chin et al. ( 1993), we 
also made use of some data which had not yet been pub- 
lished. These consisted of small grid maps made in the 
HCN (7=1 ->0) and HCO + (7=1 ->0) lines in two observing 
runs with the SEST in 1997 January and July. The maps 
cover a region of roughly 1 .'5 x 1 '.5 centered on the position 
a 2 ooo=5 h 13 m 18. s 2, 6 2 ooo=-69 22'35". The grid spacing was 
20" (for comparison the SEST beam has a FWHM of 58"), 
but the grids were not completely sampled in the outer parts. 
A 3mm SIS receiver was used in single-sideband mode with 
a 2000 channel acousto-optical spectrometer having a 43 kHz 
channel separation. 

A dual beam switching mode, where the source is placed 
alternately in the main and reference beams, was used to re- 
move the sky background, and the spectra were placed on 
a Tmb scale using a main beam efficiency of 0.76. Typical 
single-sideband (SSB) system temperatures were ~200 K. 
The CLASS package was used to calibrate and baseline the 



spectra. For purposes of this paper we have produced high 
signal-to-noise spectra for both lines by averaging the 9 spec- 
tra within -20" < Aa < 20" and -20" < AS < 20" relative 
to the center position. The corresponding maps show only 
a moderately resolved central source and are not reproduced 
here. Further details of the SEST observations will be pre- 
sented in a forthcoming paper (M. Wang et al., in prepara- 
tion). 

2.4. Mopra data 

A 4' x 4' region centered on the position ct2ooo=5 h 13 m 18 s , 
^200o =_ 69°22'25" was mapped in the 7=1^0 transitions of 
HCO + , HCN, and CO with the 22-m ATNF Mopra telescope 
in three runs from 2004 June to September. The newly imple- 
mented On-the-Fly (OTF) mode was used, in which the tele- 
scope takes data continuously while moving across the sky. 
Spectra were spaced by 7" so that the 36" Mopra beam (33" at 
CO) would be well oversampled in the scanning direction; the 
row spacing was 10", also oversampling the beam. Each row 
was preceded by an off-source integration at a location about 
10' from the map centre. Typical SSB T* s values of 200 K 
for HCO + and HCN and 600 K for CO were obtained in clear 
conditions. The pointing was checked on the SiO maser R 
Dor prior to each map; typical corrections were less than 5". 
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FIG. 2. — (a) ATCA 1.3cm map (contours) overlaid on an R-band image from the Digitized Sky Survey. Contour levels correspond to 2n 2 mJy bm ' (n=l,2,...). 
(b) Mopra CO (1-0) map (black contours) overlaid on an H I peak intensity image (grayscale and gray contours). Contour levels are 10 K km s for the CO 
map and 40, 50, ..., 80 K for the H I map. The dashed line represents the position-velocity slice shown in Fig.|3] (c) Mopra HCO + (1-0) map. Contour levels are 
spaced by 0.6 K km s _I . (d) Mopra HCN (1-0) map. Contour levels are spaced by 0.5 K km s . In each case the solid open circle represents the ATCA pointing 
position. Offsets are relative to a 2 ooo=5 h 13 m 18 s , 5 2 ooo=-69°22'25". 



The digital correlator was configured to output 1024 chan- 
nels across 64 MHz in each of two orthogonal polarisations, 
except for the June observations where only one polarisation 
was recorded. Initial spectral processing (baselining and cal- 
ibration onto a T A * scale) was performed using the LiveData 
task in AIPS++, and the spectra were gridded into datacubes 
using the AIPS++ Gridzilla task. During the gridding, a Gaus- 
sian smoothing kernel with a FWHM half that of the beamsize 
was convolved with the data, so the effective resolutions of the 
final cubes were 40" for HCO + and HCN and 36" for CO. The 
cubes were then rescaled onto a Tmb scale using an "extended 
beam" efficiency of 77 x b=0.55 for CO and 0.65 for HCO + and 
HCN (Ladd et al. 2005). This takes into account that sources 
larger than about 80" in diameter will couple to both the main 
beam and the inner error beam of the telescope. 

Comparison of Mopra and SEST spectra taken at the cen- 
tral Mopra position, with the Mopra maps convolved to the 
SEST resolution, indicate that the Mopra line intensities are 
lower by ~30%. This is not surprising given that the extended 



beam effici ency is higher th an the main beam efficiency rj m \, 
by ^30% ( Ladd et al. 2005). While use of ?7 m b rather than -q^ 
would place the Mopra data on a more accurate temperature 
scale for compact sources, it would lead to an overestimate 
of the total map flux when assuming a Gaussian beam size 
equal to the main beam. Thus we continue to use -q^ to cali- 
brate the spectra, noting however that the quoted beam size is 
somewhat less than the effective beam size due to the presence 
of a significant inner error beam. 

3. ANALYSIS & RESULTS 

3.1. The N113 Molecular Cloud 

We begin by presenting maps and derived properties for the 
molecular cloud as a whole, based on the Mop ra da ta. We 
discuss the high-resolution ATCA 3-mm data in 33.21 

3.1.1. ATCA 1.3-cm and Mopra maps 

Figure |2ja) shows the ATCA 1.3-cm continuum image 
overlaid on a red Digitized Sky Survey image showing 
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TABLE 2 
Parameters of Gaussian fits to Mopra images 



Line 


Center RA Center DEC 


Major axis a 


Minor axis" 


Position angle 




(J2000) (J2000) 


(") 


(") 


(°) 


CO 


5 13 19.9 -69 22 43 


81 ±2 


67 ±2 


-80 


HCO+ 


5 13 19.4 -69 22 32 


89 ±3 


67 ±2 


-59 


HCN 


5 13 19.3 -69 22 33 


84 ±4 


25 ±2 


-4 



"Deconvolved using a Gaussian beam of 36" FWHM for CO and 40" for 
HCO+ and HCN. 



much of the diffuse Ha emission from the H II region. In 
addition to the strin g of three bright sources reported by 
iBrooks & Whiteoakl(ll997l) . labelled 2, 3, and 4, we have de- 
tected three additional continuum sources along the same for- 
mation, one west of the main source (Source 1) and two to the 
southeast (Sources 5 and 6). The fluxes of the six continuum 
sources, from west to east, are estimated based on Gaussian 
fits to be 6.2, 66, 18, 4. 1, 3.9, and 4.4 mJy. Measuring the flux 
by direct summation of image pixels yields fluxes consistent 
with these values within 1 mJy. Apparent sources which lie 
outside this linear structure (e.g., southwest of the image cen- 
ter) are likely to be sidelobes that were incompletely removed. 
Figures Etb)-(d) show the Mopra CO, HCO\ and HCN 
maps integrated from a velocity of 23 1 to 240 km s" 1 . Re- 
markably, at the resolution of Mopra the molecular cloud ap- 
pears to be a centrally condensed structure, in contrast to the 
filamentary structure which is characteristic of the H I clouds 
in the LMC (Kim et al. 2003). The CO map is shown overlaid 
on H I emission contours from the combined Parkes/ATCA 
The HCO + emission imaged with 



survey (iKim et a l. 2003). 
Mopra shows a similar extent to the CO, however the peak 
intensity is shifted ^10" to the north. The HCN emission ap- 
pears to be narrower in the east-west direction than the other 
tracers, as confirmed by fitting of two-dimensional Gaussians 
to the images (Table [2j. 

3.1.2. Mass and Dynamical Properties 

We can estimate the mass of the molecular cloud from 
the CO data in two ways. Adopting a conversion from K 
to Jy of 1 K = 9.6 x 10" 3 Jy arcsec" 2 , the total CO(1-0) 
line flux detected by the Mopra observations is 2.9 x 10 3 Jy 
km s" 1 over a 4' x 4' region. Using a standard (Galactic) 
CO to H 2 conversi on factor of N H2 /I C q = 2.0 x 10 20 cm" 2 
(K km s" 1 )" 1 dStrong & Mattox 1996; Dame et al. 2001), this 
translates into a total gas mass (including He) of 8.2 x 10 4 
Mq. We can compar e this with a virial mass estimate, fol- 
lowing Heikkila et al1 (ll999t) . of M vir = 151fl C o(Av) 2 , where 
Rco is the deconvolved FWHM size of the cloud in pc. For 
N113, a Gaussian fit to the integrated CO intensity gives 
R = 74" (19 pc) after deconvolution, which for Av = 5.2 
km s" 1 (the average FWHM linewidth over the inner 1' x 1') 
yields M v ; r = 7.8 x 10 4 M©. Although o nly a rough est imate, 
this suggests, as noted previously by Chin et al. (199 a), that 
a standard conversion factor yields mass estimates for LMC 
molecular clouds that are consistent with virial mass estimates 
when observed at resolution s of ~10 pc. A sim ilar result was 
obtained for SMC clouds bv lRubio eiaTI (119931) . 

Using the mass estimate above and a radius of 15 pc (esti- 
mated as 80% of the FWHM, following Heikkila et al. 1999) 
gives an overall density for the cloud of («h) ~ 200 cm -3 . 
This is significantly less than the critical density for excitation 
of HCO + and HCN (~10 5 cm" 3 ), indicating that the material 
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FIG. 3. — Position-velocity slice through the Mopra CO data cube 
(grayscale with light contours) with H I contours overlaid. The H I data 
have a resolution of 60" FWHM. The slice is centered on a2000=5:13:18, 
<52ooo =- 69:22:43 and has a position angle of 107°. Contour levels are spaced 
by 0.5 K for CO and 15 K for H I. 



in the cloud is strongly clumped. Of course, our ATCA ob- 
servations (described below) provide direct evidence of this 
clumping. 

Inspection of the CO datacube reveals a velocity gradient 
of ^5 km s" 1 from the center to the southeast extension. The 
presence of such a gradient might be interpreted in terms 
of rotation (e.g., Rosolowskv et al. 2003). However, the CO 
position-velocity diagram (Figure^ is not strongly indicative 
of rotation, but rather suggests a high velocity dispersion core 
that blends into a low velocity dispersion wing with a slightly 
different velocity centroid. We note that the wing emission 
could be composed of unresolved, clumpy structure. An over - 
lapping slice through the H I datacube of Kim et al. (2003), 
smoothed to 1' resolution, shows a similar velocity gradient 
across a filament-like structure 10' (150 pc) in length. We 
conclude that the CO velocity gradient is not due to the ro- 
tation of a kinematically decoupled structure, but instead re- 
flects the bulk motion of the gas from which the molecular 
cloud formed. 

The CO cloud appears to sit within a local minimum of H I 
peak brightness temperature [Fig. 0b)]. This does not neces- 
sarily require a a large-scale conversion of the H I into H2, as 
the H I at the location of the cloud still has a high integrated 
intensity. Rather, the lower peak Tt is mainly the result of 
a flat-topped line profile, which may indicate absorption by 
colder, optically thick H I. 

3.2. The Dense Embedded Clump 

3.2.1. ATCA 3 -mm maps 

Figures [Hb)-(f) show the integrated intensity maps for the 
four lines observed with ATCA as well as the continuum at 
86 GHz. Both the line and continuum emission appear to be 
concentrated in a single emission peak, or "clump." These 
maps were produced by directly summing the velocity chan- 
nels over the ranges given in Table [2 Due to the presence 
of large negative sidelobes in the maps, likely resulting from 
gain calibration errors and imperfect deconvolution, the inte- 
grated fluxes in Tableware given for both the inner 10" x 
10" and 20" x 20" (larger values correspond to the smaller 
box), and should be considered indicative only. Also marked 



TABLE 3 
Comparison of Mopra and ATCA fluxes 



Line Mopra (tot) Mopra (pb) ATCA (im) ATCA (uv) a ma j (uv) a m \ n (uv) 



CO 


2900 


320 


— 


— 


— — 


HCO+ 


150 


17 


8.8±0.2 


14.2±0.9 


7"7±0"8 5"4±0"8 


HCN 


73 


10 


6.9±0.2 


10.5±0.9 


6"3±2"2 4"8±2"2 


HNC 


— 


— 


2.1±0.1 


3.9±0.5 


6"8±3"0 5"2±2"4 


C 2 H 


— 


— 


1.5±0.2 


3.9±1.1 


7"0±1"7 — 


Cont. 


— 


— 


42mJy 


47±7 mJy 


2"6±0"6 — 



NOTE. — All fluxes except continuum are in Jy km s . Mopra fluxes are derived from the total image (2nd col.) and restricted to the ATCA field of view (3rd 
col.). ATCA fluxes are derived from the central 10" X 10" of the image and from a Gaussian fit to the visibility data. The Gaussian has been restricted to be 
circular for CiH and the continuum. 



on each panel of Figure [2 are the positions of the obscured 
2MASS point sourc e 05131 7 77-69 22250 and the H 2 maser 
detected by Lazend ic et alJ (|2002). Note that the infrared 
source coincides with the continuum source but is offset from 
the peak of the molecular line emission. 

3.2.2. Flux Recovered by ATCA 

Table [3] lists the total fluxes detected in the Mopra and 
ATCA maps, as well as the flux detected by Mopra with the 
ATCA primary beam response applied (i.e., the Mopra flux 
restricted to the ATCA field of view). We caution that the Mo- 
pra flux within a single primary beam may be overestimated if 
emission outside the beam is smoothed into it; it may also be 
underestimated because of the use of r]^ to calibrate the data. 
In any case, we find that the ATCA flux within the central 10" 
x 10" of the integrated HCO + and HCN maps is 50-70% of 
the Mopra flux within the ATCA primary beam. For an inde- 
pendent comparison, we can also extrapolate the ATCA flux 
to short ("single-dish") baselines by fitting a two-dimensional 
Gaussian source in the visibility domain. As shown in the 
last three columns of Table [3] which give the fitted flux and 
FWHM major and minor axes, this yields a total flux that is 
consistent with the single-dish flux within the ATCA primary 
beam (assuming ^20% calibration uncertainties for both tele- 
scopes). The difference compared to measuring the flux in 
the image plane is likely due to a more inaccurate extrapo- 
lation to short baselines provided by CLEAN, coupled with 
uncorrected phase errors that lead to negative regions in the 
deconvolved images. 

The analysis above suggests that the more extended com- 
ponent of HCO + or HCN emission which does not appear in 
the ATCA images can be inferred by inward extrapolation of 
the visibility amplitude. Assuming all of the flux is recov- 
ered by the visibility plane fitting, we derive an approximate 
beam filling fraction (within the 36" ATCA beam) of 3% for 
the HCO + and HCN emission. Of course, this applies only 
to the region observed with the ATCA, which contains just 
~12% of the total flux detected by Mopra. Clearly it would 
be of great interest to cover the entire molecular cloud with 
an ATCA mosaic, to determine what fraction of the total line 
emission the interferometer would be sensitive to. 

3.2.3. Fitting of ATCA and SEST spectra 

While the Mopra OTF maps of HCN and HCO + cover a 
larger area, their sensitivity at any given position is a factor of 
3— A worse than the pointed SEST spectra, so we have opted 
to derive line widths and hyperfine ratios from the SEST data 
for comparison with the ATCA data. 
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FIG. 4. — Comparison of ATCA and SEST spectra. SEST spectra, in 7" m b 
units (top two panels), have been averaged over a 40" X 40" region centred 
on a2000=5:13:18, 52ooo= _ 69:22:35. Spectra through the ATCA datacubes, 
smoothed to 2 km s -1 and converted to brightness temperature units using 
the synthesized beam size, have been taken at the nominal peak position of 
"2000=5: 13: 17.2, 5 2 ooo=-69:22:23. Gaussian fits to the HCO+ and HCN lines 
(the latter taking into account the spacing of the three hyperfine components) 
are shown as dotted lines. 



The averaged SEST spectra are shown in Figure |4] along 
with ATCA spectra taken through the central 1" pixel of 
each datacube. The HCO + and HCN brightness tempera- 
tures Tb measured in the ATCA maps are a factor of ~10 
higher than the SEST measurements, which is not surprising 
given the presence of clumping and the factor of ^100 dif- 
ference in beam area. We have also made single-component 
Gaussian fits to the HCO + spectra and three-component fits 
to the HCN spectra; the latter assume that the two satellite 
hyperfine transitions (,F=0-1 and F=\-\) are separated by 
-7.07 and 4.85 km s" 1 from the main component and have the 
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same dispersion as the main component. We find for HCO + 
a FWHM linewidth (corrected for instrumental broadening) 
of 5.8 km s" 1 for SEST and 5.5 km s" 1 for ATCA. For HCN 
the linewidths are 4.6 km s" 1 and 5.2 km s _1 for SEST and 
ATCA respectively, though these are much more uncertain, 
being based on a multi-component fit. We also find a veloc- 
ity offset of -2.0 km s _I in the ATCA HCN spectrum relative 
to SEST, which we attribute to limited spectral resolution and 
signal-to-noise in the ATCA data. 

The ratios between the hyperfine components of HCN can 
be used as a diagnostic of the optical depth, if LTE excitation 
is assumed (e.g., Fig. 2 in Harm 1989). In the optically thin 
case, the ratio of the F=0-l/F=2-l transitions, referred to as 
R02, should be 0.2, whereas the ratio of the F=l— l/F=2— 1 
transitions, referred to as R12, should be 0.6. Fitting the SEST 
spectrum yields /?02=0.26±.01 and /?i2=0.57±.03, where the 
errors are only the formal errors on the fit and thus a lower 
estimate. These appear consistent with optically thin emis- 
sion. On the other hand, for the ATCA spectrum we find 
/?02=0.40±.12 and 7?i2=0.79±.19, indicating higher opacities 
for the dense clump at the lcr-2er significance level (note 
however that both ratios are greater than their LTE values). 
While still requiring confirmation with higher spectral reso- 
lution data and/or observations of rare isotopomers, this sug- 
gests the presence of an optically thinner HCN component 
which dominates the total flux. We suggest that, rather than 
constituting a truly diffuse interclump medium, this compo- 
nent is an enshrouding "clump envelope" that accounts for the 
difference between the fluxes measured in the ATCA images 
and those inferred by visibility fitting (Table[3j. 

3.2.4. Clump size and line width 

The visibility fits given in Table [3] indicate that the HCO + 
emission detected by the ATCA comes from a region 8" x 
5" (2.0 x 1.3 pc) in size (FWHM), with the HCN emis- 
sion coming from a somewhat smaller region. This im- 



TABLE 4 

Integrated line intensity ratios towards the N113 peak 



Ratio 



SEST a ATCA 30m ATCA 45m ATCA 75m 



HCO+/HCN 1.35±0.06 


1.27±0.09 


0.93±0.13 0.90±0.29 


HCN/HNC 2.81±0.21 


2.68±0.24 


3.49±0.71 1.67±0.47 


C 2 H/HCN 0.29±0.02 


0.30±0.05 


0.32±0.08 <0.56 



"From Chin et al. 1 1997). 



plies a radius of ^1.6 pc, comparable to that of the largest 
clumps which are found in Galactic massive star-forming 
regions (Mookerjeaet al. 2004; Reid & Wilson 2005]), and 
much larger than the typical size (~0.1 pc) of a core (a con- 
densation which forms only a single star or a multiple star 
system). We therefore believe it is appropriate to refer to this 
object as a "clum p." We n ote, however, that it is common in 
the literature (e.g., Plume et al. 1997) to use the term "massive 
cores" to refer to the birthplaces of the young stellar clusters 
where massive stars typically reside. 

It is difficult to estimate the mass of the clump from the ob- 
served line fluxes given considerable uncertainties in the ex- 
citation, abundances, and optical depths of HCO + a nd HCN. 
However, we can derive a virial mass as in 33.1.21 this time 
using R=\.6 pc and Av=5.5 kms" 1 : M vir = 151fl(Av) 2 = 
7.3 x 10 3 Mq, or about 10% of the total cloud mass infer red 
from the CO data. Alternatively. iGao & SolomonH2004al) ar- 
gue on the basis of large velocity gradient (LVG) calculations 
that the mass of dense gas traced by HCN satisfies 



M dense (H 2 ) < 25 



^HCN 



K km s ! pc 2 



Mr. 



They obtain this result assuming an HCN abundance of 2 x 
10~ 8 , a kinetic temperature of 20-50 K and an intrinsic HCN 
brightness temperature of Tj, > 5 K. Using this formula gives 
an upper limit of 1.3 x 10 4 M Q for the N113 clump, within 
a factor of 2 of the virial estimate. The agreement between 
the two mass estimates, along with the high densities implied 
by the detection of HCN, HCO + , and HNC, suggests that the 
assumption of virialization is reasonable. The average density 
of the clump is inferred to be («h) ~ 10 5 cm" 3 , close to the 
critical densities for HCN and HCO + . 

3.2.5. Line ratio analysis 

Figure [5] shows the integrated flux (in Jy km s" 1 ) for all 
four lines observed with ATCA as a function of the projected 
baseline. The equivalent curve for the continuum flux (in Jy, 
scaled down by 100) is also shown. Since the ATCA spacings 
are roughly in multiples of 15m, data have been averaged in 
bins of that width. As expected for a resolved source, the 
fluxes decrease on longer baselines, except for C2H, which 
shows roughly equal fluxes on 30m and 45m baselines (it is 
not detected on longer baselines, however). The relatively 
faster decline seen for HCO + compared to HCN is responsi- 
ble for the larger source size inferred from the visibility fit- 
ting (Table |5|l and is corroborated by the Mopra maps. For 
the shorter ATCA baselines, we have tabulated several in- 
tegrated line ratios, 7(HCO + )//(HCN), 7(HCN)//(HNC), and 
7(C 2 H)//(HCN), in Table H These are give n alongside the 
SEST line ratios reported by Chin et alJ < 19971) . Note that for 
C 2 H we only include the 87.3169 GHz line in these calcula- 
tions. 

4. DISCUSSION 



4.1. Comparison with Galactic Clumps 

Recently a number of papers have analyzed the sub- 
structure of molecular clouds hosting massive star formation 
in our Galaxy (e.g., Shirley et al. 2003; Faundez et al. 2004). 
The clumps found in these studies have typical masses of 
100-1000 Mq, sizes of 0.25-0.5 pc, and peak densities up 
to 10 6 cm" 3 . The clump we have imaged in Nl 13, with a ra- 
dius R ~ 1.6 pc and mass M ~ 10 4 Mq, lies at the upper end 
of this distribution in terms of size and mass. We believe our 
size estimate to be realistic: while the resolution of the ATCA 
maps (^4") is only slightly smaller than the deconvolved size 
of the clump (6"), the visibility amplitude as a function of 
baseline indicates that the source is well-resolved (Fig.[5j. It 
is tempting to speculate that such massive clumps may be re- 
sponsible for the formation of " populous" clusters of 10 3 — 10 4 
stars in the LMC (van den Bergh 1991 ). 

Consistent with some previous studies of massive star- 
forming regions (e.g. Plume et al. 1997), we do not find ev- 
idence of a size-linewidth relation of the form Av oc R - 5 as 
seen in low-mass star-forming regions (Larson 1981). Indeed, 
the velocity dispersion of the clump is comparable to or even 
larger than that of the molecular cloud as a whole (~5 km s _1 ). 
This may indicate a high level of turbulence, as found in mas- 
sive star-forming cores in the Galaxy. However, we caution 
that the observed line width could be affected by optical depth 
effects, and confirmation of this result with an optically thin 
tracer is required. 

4.2. The Star-forming Enviroment in N113 

The overall picture that emerges of the interstellar medium 
near N 1 1 3 is of an environment which has been strongly af- 
fected by recent star formation. In addition to the diffuse Ha 
seen in optical images, the linear sequence of 1-cm contin- 
uum sources, including two which are associated with H2O 
maser emission, point to the presence of young ionizing stars 
[Fig- Et a )]- The brightest H2O maser is found near a rela- 
tively weak (^4 mJy) continuum source (No. 4) southeast of 
the ATCA field, where there is still substantial HCO + and CO 
emission but no appreciable HCN. Moreover, referring to the 
2MASS Point Source Catalog (Skrutskie et al. 2006), the as- 
sociated near-infrared source is apparent in all three 2MASS 
bands (mj = 13.87, mn = 13.46, Mk = 12.71), suggesting it 
is relatively unobscured. These two lines of evidence point 
to an object that is somewhat more evolved than the bright- 
est (^66 mJy) of the continuum sources, which is near the 
peak of the molecular line emission and detected in Zf-band 
but not /-band by 2MASS. This obscured source (ray > 14.97, 
niH = 14.92, mK = 13.25) is shown by a solid square in Fig.^ 
A third line of evidence is the presence of an OH maser in the 
former object, again suggesting a later stage of evolution. 

Assuming the brightest continuum source emits optically 
thin free-free emission (spectral index -0.1), the measured 
fluxes at 24 and 86 GHz (66 and 42 mJy respectively) cor- 
respond to an ionizing photon output of 1.5-3 xlO 49 s _1 
(Condon 1992), equivalent to one or two 06 stars (e.g., 
iTurner et al.lll998l) . Although the measured fluxes suggest a 
spectral index of -0.3 instead of -0.1, the larger calibration 
uncertainties of the 86 GHz data (a result of higher T* and 
relying on an SiO maser for gain calibration), possibly com- 
bined with resolving out some extended flux, could explain 
this discrepancy. The 86 GHz image also appears to show 
an elongation to the north, whose significance is unclear but 
is interesting in light of the presence of an H2O maser ~4" 



north of the Tf-band source. The maser and associated contin- 
uum emission could signify an outflow from the central ioniz- 
ing source, but given their large projected distance from that 
source (^1 pc), are more likely to trace other ionizing stars in 
the same association. 

Regardless of the exact nature of the radio continuum emis- 
sion, its location only 3."6 arcsec (0.9 pc) from the dense 
molecular clump will surely affect the lifetime of the clump. 
One possibility is that the pressure of the expanding H II re- 
gion may compress the clump and thereby incite gravitational 
collapse. Whitworth et al. (1994) discuss how shocked gas 
layers generated by expanding H II regions tend to fragment 
into ~7 Mq condensations which can become gravitation- 
ally unstable, yielding a self-propagating mode of star for- 
mation. Such a scenario might explain the alignment of star- 
forming regions seen in the continuum. On the other hand, 
Gorti & Hollenbach (2002) discuss how FUV radiation sur- 
rounding an H II region heats the clumpy molecular gas and 
eventually photoevaporates it, limiting the efficiency of star 
formation. A clump as massive as that inferred by our ATCA 
observations would not be easily photoevaporated, but might 
initially experience a shock compression due to heating and 
expansion of its outer layers. Observations of shock tracers 
in this region with Spitzer may be able to directly trace this 
process and assess its ability to trigger new star formation. 

4.3. Evidence for PDR Chemistry 

The large HCO + /HCN ratio observed with the shortest 
ATCA baseline (and in the Mopra data), decreasing on longer 
interferometer baselines, is consistent with the HCO + emis- 
sion originating from a more extended region than HCN . Sim- 
ilar co nclusions for the LMC were reached by Johansson et al. 
(1994 ), based on SE ST mapping of the N159W cloud, and 
Heikkila et al. ( 1999), based o n small maps of sever al addi- 
tional clouds. As pointed out by Heikkila et al. ( 1999), two of 
the major chemical pathways for the formation of HCO + in- 
volve C + , which is expected to be abundant near the surfaces 
of clouds exposed to FUV radiation, i.e. in photon-dominated 
regions (PDRs). Thus, given the higher [C II]/CO flux ra- 
tios in the LMC (Mochizuki et al. 1994), and especially to- 
wards bright H II regions (Israel et al. 1996), HCO + emission 
likely traces lower density gas in the LMC in addition to well- 
shielded dense cores. The situation is less clear for HCN: 
while it can be produced by endothermic reactions between 
CN and H2 in the outer layers of a PDR, it is easily be dissoci- 
ated back to form CN in the presence of strong FUV ra diation 
( Sternberg & Dalgarno 1995). We discussed above ( ^3.2.3> 
the possibility that the more optically thin HCN arises from 
a "clump envelope" which is not fully imaged in our ATCA 
data. The larger size of the HCO + emission region then im- 
plies a somewhat larger envelope in this line than in HCN. 

C2H is also expected to be enhanced when C + is, and is 
especially strong relative to CO in the actively star-forming 
cloud N159W (Heikkila et al. 1999). However, the SEST and 
ATCA observations show a roughly constant 7(C2H)//(HCN) 
from large to small scales, suggesting that C2H may trace 
dense gas rather than the PDR region. Such a result is rather 
surprising, given the relatively low dipole moment of C2H, 
which should lead to its excitation in even low-density re- 
gions. An alternative explanation is that C2H appears less 
extended than HCO + because of lower optical depth. The 
hyperfine structure of C2H (N=l-0, 7=3/2-1/2) consists of 
two lines (7^=2-1 and F=l-0) separated by 12 MHz. The 
integrated intensity ratio of the 2-1 to 1-0 line measured by 
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FIG. 6. — (a) 8/xm MSX image (heavy contours) overlaid on the Mopra CO map from Fig.|5](light contours and grayscale). Contour levels for the MSX image 
are spaced by 1.5 X 1(T 10 W cm -2 sr _I , (b) 21/im MSX image (heavy contours) overlaid on the 1cm continuum map (light contours and grayscale). Contour 



levels for MSX image are spaced by 4 X 10 l0 W cm 2 sr ', while the continuum image has contours of 2ir mjy bm 
is shown as a solid circle. The position of the bright H9O maser is shown by a cross. 



(n=l,2,...). The ATCA pointing position 



SES T is 1.8 JChin et al.lll997l) . close to the theoretical value 
of 2 llTucker et alJ 19741) . and indicating that the opacity is rel- 
atively low. On the other hand, the HCO + /H 13 C O + intensity 
ratio towards N113 as measured by SEST is 34 (Chin et al. 
1996), consistent with the isotopic 12 C/ 13 C ratio derived by 
Johansson et al. ( 1994) for N159. Thus, the optical depth in 
HCO + is likely to be small as well. More sensitive observa- 
tions will be needed to draw firmer conclusions on how C2H 
is distributed relative to HCN and HCO + . 

A more direct indicator of PDR conditions may be the in- 
frared emission at 8/im, which is primarily attributed to large 
organic molecules known as polycyclic aromatic hydrocar- 
bons (PAHs) heated by far UV radiation in the outer layers of 
the cloud. As Figure|6ja) indicates, the distribution of 6.8/im- 
10.8 fim emission near Nl 13, as measured by the MSX satel- 
lite JPrice et al]ll998l) . closely matches the CO distribution 
in extent. This provides further evidence that the extended 
HCO + emission, which follows the CO distribution, is en- 
hanced due to PDR chemistry. The 21/^m emission, which 
traces heated dust, is peaked slightly to the southeast of the 
PAH emission, and corresponds well to the 1.3-cm radio con- 
tinuum emission [Figure|6tb)]. 

4.4. Comparison of HCN and HNC 

Aside from a highly uncertain measurement for the 45m 
baseline, the HCN/HNC intensity ratio appears to decrease on 
longer baselines, which sample smaller scale structure. This 
would a confirm a prediction by Chin et al. ( 1997), that the 
ratio would be large in warm gas subject to strong UV heat- 
ing but approach unity in cloud cores. The difference may be 
due to additional neutral-neutral pathways for HCN formation 
that occur in a partially ionized medium, such as CH2 + N — > 
HCN + H (Goldsmith et al.ll981l) . These produce HCN faster 
than the equivalent reactions (NH„ + C — > HNC + H n _i) pro- 
duce HNC, due to the higher proportion of CH„ radicals and 
the tendency for C to be ionized in PDRs. However, the abil- 



ity of such reactions to produce HCN/HNC abundance ratios 
much different from unity has been questioned (Herbst et al. 
2000). In any case, the observed ratios are consistent with 
the HCN/HNC ratios of ~l-6 observed in starburst galaxies 
( Aalto et al. 2002). A major uncertainty in interpreting inten- 
sity ratios as abundance ratios is the assumption of low optical 
depth. The decreasing HCN/HNC on small scales coincides 
with an apparent increase in the HCN opacity, which could 
affect the line profiles through saturation or self-absorption. 

5. SUMMARY & CONCLUSIONS 

We have presented high-resolution imaging of HCO + , 
HCN, HNC, and C 2 H towards the N113 H II region in the 
LMC, performed with the ATCA interferometer, as well as 
single-dish mapping of the CO, HCO + , and HCN lines with 
Mopra. Our main conclusions are summarized below. 

1. On scales of several arcminutes, the molecular cloud 
(M ~ 10 5 M©) shows a weak velocity gradient which is 
shared with an H I filament which encompasses it. Sev- 
eral radio continuum sources lie along the center of the 
filament, making it the locus of recent star formation. 

2. Near the brightest continuum source, which we inter- 
pret as a compact H II region, we detect a massive 
clump of dense molecular gas whose peak is slightly 
offset from the 3mm continuum peak. The clump has a 
radius of ~1 pc, and if in virial equilibrium, a mass of 
~10 4 Mq, about 10% of the cloud's mass. This would 
imply a density of ~10 5 cm" 3 . With the ionizing flux 
of 1-2 O6 stars, the compact H II region may possibly 
trigger the formation of new stars in the clump. 

3. While the clump — when combined with a more dif- 
fuse envelope inferred from visibility fitting — may ac- 
count for most of the HCO + and HCN emission within 
the ATCA primary beam, the existence of additional 
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clumps within the cloud is not excluded given the small 
extent of the region mapped thus far with the ATCA. 

4. A decreasing HCO + /HCN ratio on longer interferome- 
ter baselines, which probe smaller scales, suggests that 
the HCO + emanates in part from a more diffuse region 
then the HCN. This is consistent with the larger spatial 
extent of the HCO + emission in the Mopra map and its 
agreement with the CO and 8/im morphology. 

5. We also see evidence for a decrease in the HCN/HNC 
ratio on smaller scales. This may reflect a lower ratio 
in cold, dense regions, as found in the Galaxy, but the 
interpretation may be clouded by an increased opacity 
in the HCN line. 

These observations demonstrate that maps of LMC molec- 
ular clouds on scales of a few pc can be used to investigate 
the internal structure and chemistry of clouds for compari- 
son with Galactic samples. With the ATCA now covering a 
frequency range of 85-105 GHz with five antennas, further 
observations will soon achieve better sensitivity for a larger 
number of molecular lines, allowing for a more detailed com- 
parison with PDR models. Ultimately, the ability to observe 
both low and high excitation lines at high resolution with the 
Atacama Large Millimeter Array (ALMA) will facilitate a 



comprehensive understanding of the properties of molecular 
clouds in an external galaxy. 
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